makecov <- function(rho,n) {
   m <- matrix(nrow=n,ncol=n)
   m <- ifelse(row(m) == col(m),1,rho)  _label~rcm@
   return(m)
}
